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Abstract 

The excitation of a slow extraordinary wave in a overdense plasma from an ordinary 
wave impinging on the critical layer in the plane spanned by the density gradient 
and magnetic field is solved analytically by formulating the problem in terms of a 
parabolic cylinder equation. A formula for the angular dependence of the transmis- 
sion coefficient is derived. 
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The mode conversion of an ordinary (O) into a slow extraordinary (SX) 
wave, first observed in ionospheric plasmas [I], has lately received great atten- 
tion as part of the OXB process [2], where B denotes the electron Bernstein 
mode. This has been used with success for heating [3J, temperature diagnostic 
jl] and current drive [5] of the W7-AS stellarator and other hot overdense 
magnetically confined laboratory plasmas [6117] . 

The O-SX mode conversion has been solved in the past under the Wentzel- 
Kramers-Brillouin (WKB) approximation in one dimension (ID) [1.2.8 9|TXJ] . 
in ID in presence of a sheared magnetic field [TT], in two dimensions (2D) 
[T2"|13|ll4j and in 2D in presence of a sheared field [15J. Later it was the object 
of a numerical study [T6] . 

The present Letter provides the first analytical, non- VFKB-approximated 
solution of the problem. At the cost of a simpler geometry compared with 
Refs. [TTf 12f 13U14f 15] . it leads to a relatively simple formula which allows to 
visualize the spatial dependence (Figs. 4 and 6) and time evolution (Fig. 5) of 
the wave as it tunnels through the O-SX evanescent layer. This might replace 
time-consuming numerical models in the interpretation of direct, space- and 
time-resolved measurements of the electric and magnetic field in the mode- 
conversion region |17j . 
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Consider a plasma magnetized along z and plane-stratified along x. Suppose 
the inhomogeneity to be weak enough to use the local dielectric tensor. 

From the Vlasov-Maxwell system it follows, in the zero temperature limit, 
that 



V x (V x E) 
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where e is the dielectric tensor for a cold magnetized plasma [18 
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X = uip e /uj 2 and Y = co ce /u are the dimensionless density and magnetic 
field, u pe and u ce the plasma and electron cyclotron angular frequency. We 
seek solutions of wave equation [T] in the form: 



E(x,y,z,u) = E (x,lv) e 



ik z z 



(3) 



where the wavenumber k z is real-valued. Note that this is not an eikonal 
ansatz and that no assumption is made on the dependence of Eo on x. Note 
also, however, that k y =0, i.e. that the wave is assumed to propagate in the 
plane spanned by the density gradient Vn (direction x) and by the magnetic 
field (parallel to z). 

With this substitution, Eqjl] rewrites: 
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where N = ck/w is the refractive index. In the remainder, x will be renor- 
malized to the vacuum wavelength by replacing 
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Circular components E ± = (E 0x ± iE 0y )/\/2 permit to diagonalize the 
tensor on the right hand side: 
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The circularly polarized X-modes and the linearly, parallel to z polarized O- 
mode are indeed the independent eigensolutions of the homogeneous problem. 
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The derivatives on the left hand side introduce a medium inhomogeneity and, 
being off-diagonal, they couple the otherwise independent modes, accounting 
for their conversions. 

The diagonal elements of Eq|6] have been expressed in terms of the cutoff 
densities 

X R = (1-N 2 Z )(1-Y) (7a) 
X L = {l-Nl){l + Y) (7b) 
X P = 1, (7c) 

These correspond to different parts of the cold dielectric tensor (R, L or P, 
in Stix notation [18]) being set to 0. 

Xp is the cut-off for the "fast" (in the sense of the phase velocity) right- 
handed (R) X-mode. X L and X P , instead (Eqsj7jD-c), are the slow X- and 
O-mode cutoff densities, X S x and X Q respectively: 

X S x = X L (8a) 

X = X P (8b) 
if N z < N z>opt . If instead, N z > N z ^ opt , 

X = X L (9a) 

X S x = X P (9b) 

Here 

N l P t = p^j, (10) 

is an optimal value of N z making the O- and slow X-mode degenerate. N Z)0pt 
is the optimal value of N z yielding complete OX mode conversion [Tf2f8f9] . 
From Eqsj7 10 it can be shown that Xo is always the smallest between Xl 




and X P . 

At X ^> Xp, where no fast X-mode exists (E_=Q) we can restrict our 
attention to: 

Xl ~ X ' 1 ~f)E+ = f-^^E z (11a) 

X7i E+ (llb) 

where the subscript o has been dropped for brevity and we have changed 
coordinate from x to X. L n denotes the dimensionless local density length- 
scale subject to the same normalization as x (EqjH]). For example, L n =100 
means that it takes 100/27T vacuum wavelengths to go from the X=0 to the 
X=l location. 

Partial derivatives are now taken with respect to X, which allows not to 
make specific assumptions on the density profile X(x). The magnetic field Y, 
on the other hand, will be treated as constant because its variation in the 
density gradient region is negligible in most experiments. 
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If the cutoffs are well separated, an O-wave starts fading at the X=Xp or 
X=Xl location (whichever is smaller) and doesn't reach the other cutoff, thus 
it does not couple with the X-wave. Under these circumstances, EqsfIT] reduce 
to separate Airy equations for E + and E z , the most general solutions of which 
are linear combinations of Ai and Bi Airy functions. However, Bi solutions 
diverge and do not make physical sense here. Hence, the factor pertaining to 
the unphysical solution Bi has to be |19| . As a result, apart from factors, 



E, = Ai 



{X-X L ). 



Y + l 



(12a) 



E, = Ai 



{X - X P ) 



(12b) 



These are the solutions when the O and SX cutoff are well separated. As 
FigjTJillustrates, the O and SX waves evanesce exponentially after their cutoffs. 
Physically, this is because the polarization currents (the plasma response) 
grows up to complete cancelation of the displacement current dE/dt. As a 
result, the O-wave is negligibly small at the other cutoff and is not capable of 
exciting an X-wave. This a posteriori legitimates having neglected the coupling 
terms iN z d x /\/2 in Eqjll] 

Note in FigjTJthat the O and SX wave domains of existence are X < Xq 
and X < Xsx- This is consistent with their frequencies necessarily being 
higher than their respective cutoff frequencies, marked by arrows in Figj2^i. For 
optimal or about-optimal incidence, however, the dispersion relation modifies 
as in Figj2jD, i.e. the frequency of the mode-converted SX wave is initially 
smaller than cutoff, which implies access to densities X > X S x, up to the 
turning point [2UII212T] X = 1 + Y^f. 

It is intuitive that when cutoffs are only few decay-lengths far, an O wave 
can excite an SX wave on the other side, and viceversa. In this (general) case 
the coupling terms on the right hand sides of EqsJTT]have to be retained. 
To eliminate E z , we apply the operator (Xp — X + d\/Lfy to Eq.lli, the 



operator ijf^ to Eq. 11 d and take the difference. After some algebra, we 
obtain: 

1 1 fX L -X Xp-X iVf 
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In the vicinity of cutoffs the partial derivatives dx (proportional to iN x ) 
vanish at least linearly if not quadratically with X. After retaining only the 
zeros of lower order and taking the real part of Eqjl3] while assuming, without 
losing generality, that E z is real, we obtain: 



N'l 



l + Y 



E'! - 2E' + + L 2 n (l - X){X L - X)E + = 0. 



(14) 
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In Eqs[l3] and 14 the symbol ' denotes derivation with respect to X. In- 



stead, in terms of spatial derivatives ^ = 



dx 2 

where we have introduced 



(15) 
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are respectively an effective squared refractive index and a damping ratio 
for E+. Note that both vary with x. The spatial non- uniformity of N 2 e ^ 
implies that the wavelength changes or, when N 2 e ^ becomes negative, it 
implies that the wave becomes evanescent. The non-uniformity of £, on the 
other hand, determines where E + decreases, and how strongly. 

Unless the density gradient is extremely steep, thus L n extremely small, 
the damping ratio ( can be ignored and Eqjl5] takes the form 



d 2 E, 

dx 2 







(17) 



From Eq 



16a 



N 2 e fj can be written as aX 2 + bX + c, where a, b and c are 
constant coefficients (remember that for simplicity we are assuming uniform Y 
and slab geometry, thus, conserved N z ). Therefore, Eqjl7]can be recognized as 
a parabolic cylinder equation, the solutions of which are tabulated |22 |23(23] . 



Eq 16a can be thought of as a parabolic "potential" in X, parameterized 
in N 2 through Xl and locally representing a good approximation (Fig|3]) of 
the Appleton-Hartree dispersion relation [T8]. 
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For N 2 



;isb) 

~-N 2 opt (and only for that value) X^=Xp=l. Then the parabola in 



question develops entirely in the positive half-plane and the wave is immune 
from evanescent, N 2 e ^ <0 barriers. 

This leads to the simplest solution of Eq|l7| in terms of parabolic cylinder 
functions of order —1/2 
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[l + i)^a(X -1)]} /V2 (19b) 
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where 



2^ 
JVJ(1+Y) 



(20) 



reduces to a = 2L 2 /Y in the present, optimal incidence case (Eq, 10). The 
equivalent formulas 



C = a 1/s J\X- 1| J_i v^(X-l) 2 /2 
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(21b) 



sign(X - 1) a 1/8 * 

involve 1st kind Bessel functions of fractional order. These solutions, plotted 
in Fig|4] along with their approximations 
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aVyiX-ll 

could be termed "cosine-like" and "sine-like", for obvious reasons. Their 
linear combination according to inverse prosthaphaeresis-like formulas 



E+ = C sin uot + S cos cot, 



(23) 



with weights determined by boundary conditions at a given instant, gives 
snapshots of E + as a function of X (or equivalently, if the density profile is 
linear, of x) at different times t. Here and in Fig|5] the phase is chosen to 
yield a sine-like E + at the time origin t=0. A tendency to longer wavelength 
around X = 1 (a reminescence of reflectometric, Airy function behaviour) can 
be recognized in these graphs. Actually A — > oo in that point (and only in 
that point). This is clear because that point is a cutoff, defined by k x — > 0. 

For non-optimal incidence X^ ^ Xp and it is convenient to introduce the 
intermediate density and the semi-difference: 



X m = (X L + X P )/2 = l + 6X 



(24a) 



5X = (X L - X P )/2 = -(1 + Y) AiV 2 /2 (24b) 
where SX is proportional to the thickness of the barrier and the systematic 



N 



iV 2 t , is related to the 



error in the squared refractive index, AiV 2 
systematic error in the launching angle. 

For finite SX, the parabolic functions involved in Eqsjl9] generalize as fol- 
lows 122112311241: 



E + 
E 4 



-l + i)^a{X -1-5X) 
'l + i)^a(X -1-5X) 



(25a) 
(25b) 
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and are plotted in FigJSJ When incidence is far from optimal, Airy func- 
tions are re-obtained (Figjf^i and e). Improved launching conditions thin the 



evanescent barrier. The coupling terms in Eqs, 11 start playing a role and each 
Airy function acquires an oscillatory nature beyond its respective cutoff as a 
result of the interaction with the other (Fig|6]D and d). Optimizing N z sup- 
presses the barrier, maximizes the mutual influence between the O and SX 
solution and, thus, transmission (FigJ^). 

Transmissivity through the barrier can be defined as the ratio between the 
amplitudes at the barrier edges, i.e. at cutoffs, squared: 

This is plotted in Fig|7]and is in reasonable agreement with Mj0lus formula 

T = exp [-TrL n VZY(l + Y)(N Zt0pt - N z ) 2 ] , (27) 

here rewritten for the present geometry (N y =0) and normalizations (with 
L n corresponding to k L n of Mj0lus). In turn, this was validated against 
full wave calculations [T6f25] and agreed within error bars with experiments 
[5Rl2imtIj. In particular the two curves in Fig jT] exhibit the same width at half 
maximum. They also practically exhibit the same dependence on L n (note 
the renormalized quantity on the horizontal axis), with the angular tolerance 
augmenting if the density gradient steepens. 

In conclusion in the present Letter the problem of the mode conversion of 
an ordinary into an extraordinary mode has been formulated in terms of the 
parabolic cylinder Eqfl~7"} Their solutions for optimal (Eqs, 21 ) and non-optimal 



(Eqs 25 ) incidence on the cutoff layer allow to visualize the wave behavior in 
the vicinity of and across the degenerate cutoff (Figjl] and [5J and evanescent 
barrier (Figj6|), respectively, and with realistic amplitudes, in agreement with 



a well- validated expression for transmissivity (Fig|7]). 
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Fig. 1. O-mode (blue, solid) and X-mode (red, dashed) solutions of EqsfTT] when 
coupling is neglected, for L n =100 and Y=0.9, in the vicinity of their respective 
cutoffs. The colored stripe represents the evanescent region in between. 




Fig. 2. Cold plasma dispersion relation of ordinary (O), fast (F) and slow (S) 
extraordinary (X) waves propagating (a) perpendicularly and (b) with optimal in- 
cidence with respect to a magnetic field of strength Y=0.9. The dashed line cor- 
responds to propagation in vacuum, for reference. Arrows mark the distinct and 
degenerate, respectively, O and SX cutoffs. 
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Fig. 3. (a) Appleton-Hartree dispersion relation for optimal and slightly non-opti- 



mal propagation in Y=0.9 and (b) comparison with Eq 16a in the vicinity of cutoffs. 
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Fig. 4. Odd (purple) and even (green) exact (solid) and approximate (dashed) 
solutions of the parabolic cylinder Eq]17[ for optimal incidence. 
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Fig. 5. Time evolution of a mode-converted O-SX wave crossing the X=l cutoff 
with optimal incidence. 
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Fig. 6. Parabolic cylinder solutions, Eq 25 for values of AN^=N^ — N^ opt corre- 
sponding to (a) far too perpendicular, (b) slightly too perpendicular, (c) optimal, 
(d) slightly too shallow and (e) too shallow incidence on the evanescent barrier 
(green) . 
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Fig. 7. Transmissivity inferred from electric field values at the barrier edges in FigjH] 
(blue) and comparison with Mj0lus' transmissivity (red). The latter preserves its 
shape, the former nearly preserves it even when the density lengthscale varies by 
several orders of magnitude. 
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